Breaking Diffeomorphism Invariance and Tests for the Emergence of Gravity 
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If general relativity is an emergent phenomenon, there may be small violations of diffeomorphism 
invariance. We propose a phenomenology of perturbatively small violations of general relativity by 
the inclusion of terms which break general covariance. These can be tested by matching to the 
Parameterized Post Newtonian (PPN) formalism. The most sensitive tests involve pulsar timing 
and provide an extremely strong bound, with a dimensionless constraint of order 10 -20 relative to 
gravitational strength. 



1. INTRODUCTION 

Most approaches to new theories invoke higher gauge 
symmetries beyond those seen at low energies, with the 
known symmetries being respected at all energies. How- 
ever, it is also possible that the gauge theories in nature 
are emergent, low energy manifestations from a funda- 
mental theory with quite different degrees of freedom and 
perhaps without an intrinsic gauge invariance. Emergent 
phenomena are common in nature. While there may be 
obstacles to emergent mechanisms for fundamental the- 
ories, the possibility certainly deserves study. 

Gravitational physics is a good candidate for an emer- 
gent theory because of the poor high energy behavior of 
general relativity. While the low energy theory forms a 
good quantum effective field theory [l| , at high energies 
the perturbative theory falls apart. This may signal the 
need for new degrees of freedom and new interactions 
beyond the Planck scale. In some theories, those with a 
background independence or explicit general covariance, 
the emergent gravitational theory will fully respect diffeo- 
morphism invariance 0. However, other candidate may 
not have this feature. There are several attempts to pro- 
duce emergent gauge theories and gravity [3j , with many 
using ideas based on condensed matter emergent analo- 
gies, as well as recent work on Horava-Lifshitz theories^. 
Indeed, the Witten- Weinberg theorem Q implies that if 
emergence is to explain all the gauge theories, spacetime 
and Lorentz invariance may need to be emergent. 

Emergent theories that start from a framework without 
diffeomorphism invariance will leave behind an imprint of 
the lack of this invariance. By definition, these theories 
reduce to general relativity at low energies, plus small 
corrections. At the very least, loop diagrams probe the 
very highest energies and will be sensitive to the lack of 
the invariance at the fundamental scale. 

Gravity is also a good place for emergent phenomenol- 
ogy. The gravitational interactions are weak, suppressed 
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by powers of the Planck scale. Therefore, it is plausible 
that the residual effects of non-invariance will be rela- 
tively more visible in gravitational interactions. 

Motivated by these considerations, we will initiate a 
phenomenological study of potential small breaking of 
diffeomorphism invariance. While we do not know the 
underlying theory and therefore do not know the magni- 
tude and form of the symmetry breaking, we study an ef- 
fective Lagrangian with terms that break diffeomorphism 
invariance, and proceed to match on to the Parameter- 
ized Post-Newtonian (PPN) framework in order to pro- 
vide a remarkably stringent bound. 

The plan of this work is as follows: In Sec. 2, we 
discuss the infinite set of operators which violate diffeo- 
morphism invariance and describe a set of terms which 
involve two derivatives of the metric which we will use 
as our test cases for this study. In Sec 3, we look at the 
linearized theory, find the propagator and calculate the 
bending of light. Because gauge non-invariant theories 
have extra degrees of freedom and could have ghosts we 
use the propagator analysis to probe for the existence of 
ghosts. We find that certain values of the parameters 
are required in order to be ghost-free. In Sec. 4, we 
start the matching to the PPN framework, which is car- 
ried out in detail in Sec. 5-8. In Sec. 9 we describe the 
resulting phenomenological constraint. Sec 10 is a brief 
summary. Several appendices describe auxiliary features 
of our treatment. 



2. FORMULATION 

Once one opens up the action to include non-covariant 
terms, the possibilities are myriad. When we organize 
the theory in an energy expansion, the action is ordered 
by powers of derivatives. In a theory with a metric, the 
metric g^ix) is the primary field. There are no Lorentz 
invariant combinations of the metric without any deriva- 
tives, aside from the cosmological constant. In this case, 
the leading possible non-covariant terms in the action 
start with two derivatives. Terms with four derivatives 
or more would be suppressed at low energy. In this paper 
we study the two derivative modifications to the action. 
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One should recognize that there might be other ways 
to test diffeomorphism breaking. For example, one can 
also include terms which violate Lorentz invariance [ij], as 
such breaking is also a form of diffeomorphism violation. 
Indeed there might be good reasons for including such 
terms for an emergent theory Q. The breaking of Lorentz 
invariance has been already been studied in gravity [f| 
and we will study covariance breaking terms which are 
Lorentz invariant. 

Moreover, if one gives up covariance and gives a spe- 
cial weight to a flat metric, it is possible that one could 
consider h, v (x) — g^u{x) — rj^ as the primary held. In 
this case, there could be terms with zero derivatives in 
the actions with the leading effect being the Pauli-Fierz 
mass term Q . This term has been heavily studied and is 
phenomenologically ruled-out at essentially any magni- 
tude. When the mass is bigger than a critical fraction of 
the curvature, the van Dam-Veltman -Zakharov (vDVZ) 
discontinuity (Tol . fll| says that the predictions drasti- 
cally disagree with general relativity. When the mass is 
smaller than the critical curvature the mass generates an 
intrinsic instability 12] in the spacetime in flat, de Sitter 
or Freidman- Robertson- Walker cosmology. Only anti- 
de Sitter spaces escape these serious problems. However, 
anti-de Sitter space does not appear to be selected in 
Nature. So it appears that this form of diffeomorphism 
breaking must be identically zero. Combinations with 
higher powers of /i M „ (x) and zero derivatives should also 
be studied. Some of these are listed in Appendix A. How- 
ever, we will turn our attention to the next order in the 
derivative expansion. 

In the derivative expansion the next terms that would 
occur would be those with two derivatives of the met- 
ric. These are the ones that we study below. In addition 
there could be others with four derivatives of the metric. 
In principle these would be suppressed at low energy since 
the derivatives turn into factors of the graviton energy. 
Since we have no prior knowledge of the mass scale ap- 
pearing in the energy expansion, these should be studied 
as well, but we reserve this for future work. We provide a 
classification of the symmetry breaking operators in the 
linear approximation up to sixth order in Appendix A. 

Given these possibilities, we do not attempt a fully 
general analysis, but will look at some possibilities which 
have not been studied before and for which we can obtain 
a particularly tight bound. More general possibilities will 
be considered in future work. 

In this section, we introduce a general second deriva- 
tive Lagrangian involving the connection in ways that 
break the diffeomorphism invariance. This is a purely 
dynamical metric theory of gravity which assumes that 
|13j | that there exists a symmetric metric. In addition, 
we also assume that all non-gravitational fields couple 
universally to the gravitational field. 

Our general action takes the form 



where 



S = 



(1) 



c = 



16ttG 



(2) 



The first term above is the usual Einstein-Hilbert 
term, £ m is the matter Lagrangian, while £i are the 
diffeomorphism-violating pieces 
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(3) 



and {a,i\ are small coefficients. The Greek indices run 
over four spacetime coordinates, and we take the Lorentz 
signature to be mostly positive. Notice that £eh can be 
written as — £\ + £2 + £& — £7- Using integration by 
parts and the identities 



8 o P7 — —a pa r J - 
u vg — y 1 V a 

we obtain 

£q = £2 ~\~ surface term , 



9^,8,^9 = ^9^ (4) 



£7 — 2£i — £2 



rfacc term . (5) 



Hence, £ 6 and £7 are not independent and we drop them 
in our analysis. Moreover, using the same identities, one 
can show that any other term of square derivative, e.g. 



-gd 7 g a,3 d~ <g a p 



-gd a dpg al3 , can be uniquely ex- 



pressed as a linear combination of £\ through £5 mod- 
ulo surface terms. This set is not unique beyond lin- 
ear order, as once one has given up covariance one can 
have an infinite set of Lagrangians of the form {\f—g~) n £ 
by adding extra powers of (y/— g) n starting at the same 
two derivative order. However, such modifications do not 
contribute new terms to the linear equations of motion. 

The use of the connection in Eq. [3] is important for our 
matching to the PPN formalism. As we will see in the 
next section, this basis is larger than needed for the linear 
analysis keeping only the leading term in the expansion of 
the metric. However, the PPN formalism is sensitive to 
the higher order terms in the expansion of the metric. By 
using the connection in our operator basis, these higher 
order terms are part of expansion of the connection. 



3. LINEARIZED VERSION OF THE 
EQUATIONS OF MOTION 

Writing g, v = rj^ + h^, and retaining only the 
quadratic contribution of h, we obtain 
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where 

Ti = d 1 h aP d~ 1 h a P , T 2 = d 1 h af} d h a ~' 
T 3 = d a hd h af> , T 4 = d a hd a h. 

Inspection of eq. ((6|) reveals that the quadratic La- 
grangians 4 2) to £f are independent. Hence, we con- 
clude that the Lagrangians C\ to £4 are also independent. 

Variation of the quadratic version of the Lagrangian 
(fTj) results in the linearized equations of motion 

(-1 -ai + 3a 3 ) Uh afi + (1 + 01 - a 3 + 2a 4 ) {d a d 1 h 01 
+d d 1 h a ^) + (-1 + a 2 - 2a 4 ) rfrtdpdvhr 1 
+ (-1 + a 2 - 2a 4 ) d a d fJ h + (1 — a 2 + a 4 + a 5 ) ?7 Q ^n/i 
= 167rGT Q/3 , (7) 

where h = h%, and T a " is the energy momentum tensor 
that results by varying the matter action with respect to 

9a/3- 

Throughout this work we assume that the energy mo- 
mentum tensor of the matter fields is conserved. This 
translates in the linear theory into d a T a p — 0. Taking 
the derivative of eq. (J7J) we obtain 

2 (03 + 04) ad p h al3 + ( ai +a 2 - 03) d a d li d v h^ v 
+ (a 5 -a±)d a Uh = 0. (8) 

This equation is satisfied provided that we have either 

• (i) Trivial case: all the coefficients a 4 to 05 are set 
equal to zero, or they satisfy a 4 = 05, 03 + a 4 = 
and ai + a 2 — a 3 — ( the last two conditions 
result in ai + a 2 — 3a 3 — 2a 4 = 0, see the com- 
ments below eq. (f!5]l). This is the case of a 
diffcomorphism-invariant theory (in this case gen- 
eral relativity GR). 

• (ii) General case: all coefficients ai to 05 are dif- 
ferent from zero. In this case the realization of eq. 
(|5]) can be guaranteed if we impose the constraints 

d a h af3 = 0, and a 4 = a 5 . (9) 

The physical consequences to linear order (light 
bending) of this case is worked out below. 

• (in) Special case I: all coefficients are set to zero 
except a 4 =/= 0. In this case we make use of the field 
redefinition h af) = h af) - -q al3 h/2 to bring eq. © 
to the form 

-Uh at3 + (1 + 2a 4 ) (d a d J h l3 ' y + d^d^') 

-(1 + 2a i )r l aP d^d u h^ = 167rGT af3 . (10) 

Then, imposing the condition d a T a ^ = we get 
the constraint d a h al3 — 0. Using this constraint in 
eq. (Unj we obtain Dh a ' 3 = 16ttGT q/3 . Hence, to 
linear order, setting a 4 ^ does not lead to any 
physical consequences beyond GR. 



• (iv) Special case II: all the coefficients are set to 
zero except 05 ^ 0. In this case one can use the 
transformation h a @ = h a/3 ~ ri al3 h/A to eliminate 
05. Again, this leads to no physical consequences 
on the linear level beyond GR. 

It is important to notice here what has been stated in 
the literature that if we want a more general equation 
than Einstein's equation that respects Lorentz symme- 
try and reduces in the weak field limit to a second order 
equation, then we have to include other elements that are 
unrelated to the metric tensor or its derivatives, and we 
must give up the possibility of deriving Newton's theory 
as a limiting case (see e.g. The above statement 

is true only if we do not give up the diffeomorphism in- 
variance as a fundamental symmetry of the underlying 
manifold. However, we have shown that breaking this 
symmetry can still result in a second order equation that 
has Newton's theory as a limiting case (see the discussion 
below) provided that we take the covariance breaking co- 
efficients to be small enough. 

Also, one may argue that we should run into trou- 
bles once we break the diffeomorphism invariance. The 
classical example is Pauli-Fierz massive gravity 0. The 
argument is that since diffeomorphism symmetry is a 
dynamical symmetry, breaking it will excite the scalar 
modes that become strongly coupled even when we send 
the graviton mass to zero. This is the famous van Dam- 
Veltman-Zakharov (vDVZ) discontinuity [l(J EH- As we 
show below, this kind of discontinuity does not happen 
in our case. A simple explanation can be given here (see 
e.g. the first reference in [12] )■ The equation of motion 
for Pauli-Fierz massive gravity reads 

-ah af3 + {d a d^ + d p d 1 h a ' r ) - n a0 d»d v h fl ' / 

-d a d p h + r) al3 nh + m 2 (h al3 - r) a0 h) 

= 16ttGT q/3 . (11) 

Taking the divergence and trace of the above equation 
results in the five constraints 

d a h^ = d^h,h=^-T. (12) 

Along with the obvious behavior of 1 /m 2 in the limit of 
zero mass, the first constraint in (TT2"1) does not reduce 
to a guage condition consistent with equations of motion 
in the limit m = 0. This behavior, which signals the 
presence of a problem, is absent in our case since the 
constraint d a h a/3 — reduces to a gauge condition as we 
send Oij a 2 , etc. to zero. 



3.1. The propagator and bending of light 

To further study the linearized theory, it is instructive 
to write down the graviton propagator. This can be ac- 
complished by writing the quadratic Lagrangian in the 
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form /i M „0 MI/,a/3 h a p, and finding the inverse of the oper- 
ator O. To impose the constraint d a h a p = 0, we insert 
in the Lagrangian a term A {d a h a f}) 2 , and take the limit 
A — > oo at the end of calculations. Hence, we find 

Dp V . pcr {k) = -Arj^ripa/k 2 + B (j] w r\ V(y + r) fia i] l/p ) /k 2 
~\~A (j]p.vkpk a + Tjp a kpk v ) I fc 4 

B {j]ppk v k a + T\p a k v kp + Tj va kpkp 
~\~Tji/pk^ka) /k 4~ Ckpkvkpka / k® , (13) 

where the constants A, B, C, and D are given by 

1 — a 2 + 2<i4 



.4 
B 
C 



(1 + ax — 3*23) (2 — ai — 3a2 + 3a 3 + 60,4) 
1 



2(l + ai - 3a 3 ) 

1 — ax — 2a2 + 3a3 + 4a4 
(1 + ax ~ 3a 3 ) (2 — 04 — 3a 2 + 3a 3 4- 604) 



(14) 



and we have put 04 = 05. As expected, the propagator is 
continuous to the GR result as we set a\ = 0,2 = ■ ■■ai = 0. 

Now consider two particles with conserved energy mo- 
mentum tensor Tff. and interacting via the exchange 
of a graviton. The scattering amplitude is given by 



GT^(k)D l/tVjPtT (k)T^(k) 
G ( 

= ^2 {- AT W T (2) + ZBT^T^pv 



(15) 



The scattering amplitude between to chunks of non- 
relativistic matter is proportional to 



(16) 



where G cS = 2G(2B - A). Now, taking T?X and 
to be respectively the energy momentum tensor of the 
sun and photon we obtain the scattering amplitude, and 
hence the light bending effect 



G e ff 2 — a\ — 3a2 4- 3a 3 4- 604 
2k 2 1 — ax — 2a2 4- 3a 3 4- 4a4 



(17) 



In appendix D, we obtain the same results using the PPN 
formalism. 



3.2. Ghost Analysis 

Finally, we consider the issue of a possible ghost in- 
stability in diffeomorphism-violating theories of gravity. 
Generally, ghosts can appear in effective field theories. 
However, they usually show up only at energies above 
the cut-off scale of the theory. A thorough analysis of 
the linearized version of the action {T]) to search for such 
instabilities was given in [l5j in the context of transverse- 
diffeomorphism theories of gravity fl6| . Still, one can 
learn about these instabilities by studying the behavior 



of the momentum space propagator. A ghost propagator 
will have the opposite sign of a healthy degree of free- 
dom. Hence, existence of an instability can be read from 
(|13[) as terms with the wrong sign. A quick way of doing 
this is by sandwiching the propagator between the energy 
momentum tensors of two sources as we did above, and 
then projecting out the transverse traceless part of the 
spin-2 particle. Hence, the scattering amplitude can be 



written as 



2GB 
GB 



J (1) J (2)A^ 

a\ + &2 - 



1 

~ 2 
3a 3 



T (1) T (2) 

2a 4 



, . .) ., ■> ,. -T{i)T{2) ■ (18) 

k l 2 — ax — 3a 2 4- 3a 3 4- 6a 4 w w 

The first term is the usual spin-2 graviton coupled to 
matter, while the second term represents massless inter- 
action between conserved sources. This massless degree 
of freedom has a healthy kinetic term (not a ghost) pro- 
vided that ax + <22 — 3a 3 — 2*24 < 0. The saturation of 
this inequality decouples the massless mode. 



On the other hand, the analysis drawn in 15j showed 
that the Minkowskian vacuum in the linear version of ([T]) 
admits a linear classical instability for the vector modes 
unless we impose the constraint a 3 4- 04 = 0. These 
modes couple to the derivative of the energy momentum 
tensor, and hence they do not show up in (| 18[) since we are 
considering conserved sources. The constraint a 3 + 04 = 
along with the saturation of the above inequality (ai 4- 
a 2 — 3a 3 — 2o,4 = 0) restores the diffeomorphism invariance 
of the theory (GR) 

We show below that experimental tests of gravity 
bound the values of {a.;} to be of order 1CP 20 . This 
suppresses any significance for ghosts or other classical 
instability, if any, at low energies. At higher energies, 
one expects higher order terms, possibly higher deriva- 
tive Lagrangians, to contribute new degrees of freedom 
rendering the UV theory ghost free. 



4. NONLINEAR EQUATIONS OF MOTION 

The linear approximation is not sufficient for the PPN 
analysis that we are about to undertake. The iteration 
that is involved in that formalism mixes different powers 
of the linear field hp V . Because our original Lagrangian 
was written in terms of the connection, we are able to 
include higher order terms using the original operators. 
Moreover at this stage we are going to restrict our treat- 
ment to one of the operators in our basis, £ 3 . While 
we plan to report on it in future work, the treatment of 
the general case is very cumbersome and the use of this 
operator is sufficient to identify the strongest test, com- 
ing from the preferred frame parameter a 3 , and obtain a 
very strong constraint. 

The nonlinear equations of motion of the system can 
be found using the Euler-Lagrange equations 



0, 



(19) 
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where the comma denotes ordinary derivative. 

Taking into account the identities g a p9^ — ^2' an< ^ 
d^J—g / dg a j3 = y/—gg al3 /2 we obtain (from here on, we 
drop the subscript in 03 to reduce notational clutter) 



To obtain the Newtonian limit, we need to solve for 
500 to 0(2). Assuming g o — > far from the system we 
find goo — 2G c sU . Since we work in units in which the 
measured gravitational constant is unity, we set 



R„u - ^g^R + aM„„ = 8vrG T„ v , (20) 

where = B^+V^ , and the functions B^ v and 
are given by 



B^= - Ig^g^g^W^+g^g^g^g^T 



2^ a up 3 a " ■yc^ Sr) 

v -r A J a -a- A' 



7 A PS 

(21) 
(22) 



G e ff = 1 



(27) 



and 



The post Newtonian corrections to the propagation of 
light may be found by solving for to 0(2). To this 
order, one can use the linearized equations of motion iff]). 
However, for more involved experiments like the perihe- 
lion shift of Mercury, we need to know goo to 0(4). To 
this order we work out the full PPN parameters using 
only £3 as an example. 

In the following, we assume that the matter content is 
idealized as a perfect fluid, and hence the components of 
the energy momentum tensor to the relevant order are 



A^ = g af) g^Tl + g^g lu r 



a/3. 



(23) 



Since the equation of motion (120|) is not invariant under 
general coordinate transformations, the existence of a so- 
lution requires that we impose a consistency condition, 
as we did in the linearized case. As usual, we assume that 
the energy momentum tensor is conserved, i.e. it satisfies 
V^T^;, = 0, where V denotes the covariant derivative. In 
addition, we have from the geometry V M G M „ = 0. Hence 
we must also have 



0. 



(24) 



which is the consistency condition. 

In the rest of the paper, we use the parametrized 
post Newtonian (PPN) formalism to bound the numeri- 
cal value of a by comparing the outcomes of our theory 
of gravity with the experimental data. 



5. THE PPN FORMALISM 

To compare the different theories of gravity to exper- 
iments, we take the slow motion and weak field approx- 
imation. Such treatment is perturbative in nature and 
known as the post Newtonian formalism (l7| . In this 
method, one expands in a small expansion parameter 
which is taken to be the velocity v of a fluid element 



U 



p/ P ~n~0(2), 



(25) 



where U is the Newtonian potential, p is the pressure of 
the fluid, p is its rest mass density, and n is the specific 
energy density (ratio of energy density to rest-mass den- 
sity). The power of velocity v is 0(1), U is 0(2), n is 
0(2) and p is 0(4). In addition, since the time evolution 
of a system is governed by the motion of its constituents, 
one has d/dt ~ v ■ V, and hence 



\d/dt\ 
\d/dx\ 



0(1). 



(26) 



T 



no 



p(i + n + 

T lJ = pv l v J +p6 l 



2U) 



T w = pv 



(28) 



In addition, the metric will be constructed out of few 
gravitational potentials U, Uij, V i7 W,, $vFj 3>i> $2, $3, 
$4, $5, A, and B. The reader can refer to appendix B 
for the explicit form of these potentials as well as the 
important differential relations they satisfy. 

In the next section we review the setup used to solve 
systematically for g^ u up to 0(4). 



6. SETUP 

Far from the system under investigation, we expect 
that the metric tensor reduces to that of Minkowski 
space. Therefore, we expand the metric g^ about the 
Minkowskian background r)^ v = { — 1,1,1,1} in powers 
of v 2 



goo 



-1- 



(2) 
300 



(4) 

goo 



(3) , (5) 

goi — got + goi +■ 



9ij = 9ij + 9ij +■■■ 

(29) 



where the Latin indices run over the spatial dimensions, 

and g$ is of order v N . If we define the inverse metric 
as 



9 



= -1- 

(3) 

= 5 01 



(2) 

g Q0 

(5) 

-g° l 



(4) 

5 oo 



(2) 
n'-i 



(4) 



then using the identity g a 0g 



01 



(2) 

g 00 - 



(2) 

5oo • 



(2) 



(2) 

9ij 



52, we find 

(3) 



0i ( 3 ) 

9° =So* 



(30) 



(31) 
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The components of the affine connections are given in 
appendix C by (|C1[) . Also, the components of and 
M.P.V have the form 





(2) 


(4) 




(2) 


(4) 


-Roo 


= -Roo - 


t- Roo +■■ 


■ , Moo 


=A"loo 


+ Moo 




(3) 


(5) 




(3) 


(5) 


Roi 


= Roi -\ 


- Roi +... 


, M 0i - 


=.M« 4 






(2) 


(4) 




(2) 


(4) 




= -Roo - 


f -Ry +.. 


. , Mij : 


=Mij 4 


- Ma - 



Plugging eq. (|C2|) into eq. (|33)) and using (f36|) we find 



(1 - 3a) V 2 S= -WttG^-^ Too 



,l + 3o i 2 ) 
2 + 3a 



(2) (2) (2) (2) 

ffoo, y 4- 3fcfc,ij +(1 - 3a) g o,kk= 167rG T 

(2) _l (2) q ^ (3) is ^t?' 

ffoo, Oi 4- gkk,oi 4-(l - 3a) g i,kk= -lbirG i j . (37) 



(2) 

(32) The solution of <7qq is given by 



At this stage, it is more convenient to take the trace 
of PD)) and rewrite it in the form 



(2) 2G(l + 3a) tt 

goo = — — » — — ■ 



(l-3o)(l + 3a/2) 



(38) 



R^v 4- a FM M „ - ^g^vM J = 8tt G f T m „ - ^g M „T 

(33) 

Plugging the expansion (|32|) into eq. (j33|) we obtain, to 
the relevant order, 



Because the PPN formalism works in units in which the 
gravitational constant is unity, we must set 



G(l + 3a) 



(l-3o)(l + 3o/2) 



= 1, 



(39) 



(2) n ( (2) (2) \ ( 2 ) 

Roo +- .Moo 4- Mu J = 4^G T 00 



Hence, the normalized solutions of ggo , SiJ j an d 9oi 
takes the form 



(2) / (2) 1 (2) I (2) \ (2) 

Rij +a Mij +-6 i:i Moo Mkk = 4ttG^ T 00 



(3) (3) (3) 

R 0l +a M 0i = -8ttG T° l 

(4) a / (4) (4) (2) (2) (2) (2) " 

-Roo +2 I Moo + M u - gooMu - g, ,M, , 

( ^ (4) (2) (2) \ 

= 4vrG T 00 + I* -2 . goo T 00 , 



(34) 



where the summation is indicated in Ma- 

Finally, eqs. (|34[) have to be supplemented with the 
constraint ([24]) to every order, i.e. 



. (2) 

V 1 Mij= 

n (2) . (3) 

V° Moo 4-V 1 M 0i = 
V°7W 0l +V*M«=0. 



LOWER ORDER SOLUTIONS 



(2) 

In this section, we investigate the Newtonian, g Q0 , and 

(2) (3) 

the first post Newtonian, g\- and flfo , limits for our 
theory of gravity. By direct calculations, we obtain the 
lower order expressions for R and A4 as given by (|C2[) in 
appendix C. In addition, the first two constraints in (|3"5)) 
results in 



(2) (3) (2) 

9ij,j= ; Sto.i — 3oo,o= , 
which greatly simplifies the subsequent analysis. 



(36) 



(2) 

goo = 2U 



(2) 

9ij 

(3) 

9oi 



— 07^ + ^) 



1 



l + 3a 



[(3 + 6a)V,+ 



(40) 



8. HIGHER ORDER SOLUTIONS 

Finding goo to 0(4) is a cumbersome step since it in- 
volves all the lower order solutions . Moreover, the re- 
quirement that the solution satisfies the constraint (|2~4"|) 
makes it a long and tedious procedure. 

The solution of can be obtained from the last eq. in 
(j34|) . This involves the higher order perturbations of Roo, 
Moo and Mu. The expressions for these functions are 
given in eq. (|C3[) . Although we do not need the explicit 
value of g^j in the PPN formalism, the appearance of 
(35) this term in M^ necessitates, in general, a simultaneous 



solution for g\j and <?qq . 1 In turn, this adds more to 
the complexity of the problem by forcing us to feed the 
system in (|3~4"f with the i — j equation to 0(4). However, 
since the parameter a is small, we will be interested only 
in solutions to first order in a. To this end, one can 
solve for the terms g^' to the zeroth order of a using 

only the general relativity (GR) part, i.e. using Rf and 

neglecting completely the contribution from M^3 . At 



1 In fact, it is clear from eq. JC3I ) that we only need the combina- 
.. („ W „ W 

tion I i gu.kk ~ l 9ik,ik 
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the end, we substitute the result in (|34|) when trying to 
find <?o?- 2 This introduces an error of O (a 2 ) in our 
calculations. The solution of <?i?in GR was first given 
by Chandrasekhar and Nutku [18| in the PPN gauge. 
However, to be consistent we should get an answer that 
obeys the constraint (|36|) . We work out the details of 
these calculations in appendix E. Using the results of 
appendix E along with the dictionary of appendix F we 
find 



(4) 



(4) 



(4) 



= V 



1 (4) 1 

~2 3oo ~t; u 



7/2 + 3a 



3 a 



l + 3a 



(A 



2~ ' 1 



3a 



1 (4) 

2 3oo 



35 rr2 11 * 



2(1 + 3a) 



23 * 



16 

-2A + 2S + 



1 

Tg 



8 



6$i 



69 



$2 +8$ 3 



V 



(41) 



where V = E^, and -Ei are arbitrary functions associated 
with gjtj' as explained in appendix E. Moreover, the r.h.s 
of the last equation of (|34|) reads 



/ (4) 

AttG T 00 



(4) m (2) N 

T n -2 g 00 T 00 



= { l--a 



xV 2 [-2$i + 2$ 2 - $3 



3<J> 4 



(42) 



Now using the last eq. in (j34|) , and solving for g^ to 
the first order in a we obtain 
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gtl = I + — " 1 r-' + 



1- -a ) $w 



4 -y a 



„ 59 

3 a 

2 



+2$ 3 + (6 - 24a) $ 4 - -a (.4 + 6) 



(43) 



Next, we impose the constraint (|2~4"|) . and derive the 
equation that determines the function V. 



8.1. Determining V 

To find the condition that determines V, we use the 
constraint ([Ml) to 0(4), as in the last eq. of (|35[) . Tak- 
ing the derivative of the aforementioned equation with 



respect to x % , and using eqs. (|E1|) . (|E4|) . and (|E5[) we 
obtain 



1. 



(3) 



V'V + -V-.V,, .VI,, ,, - ( rj^w + T k kj Mij 



(2) ( 2 ) (2) (2 ) 



(2) (2) C?) (2 ) (2) (2 ) 

+ gjkM i3 . k - r£ M M - T^M ( 



(4) 



(4) 



(44) 



where S^j is given by (|E7|) , and the rest of the quantities 
are given in appendix C. Using the dictionary in JF2J, 
and solving for V we find after long, yet straightforward 
calculations 

83 rr2 5^ 45^ 7^ „^ 33^ 

V = U 2 $w $i + -$ 2 - 2$3 H $4 

16 4 w 8 2 JT 4 



(45) 



In the next section we read off the PPN parameters 
and constraint the value of a. 



9. PPN PARAMETER VALUES AND 
INTERPRETATION 

To extract the PPN parameters one has to bring the 
the form of the metric to the standard PPN metric by 
means of a gauge transformation. However, since we are 
dealing with a theory that breaks the diffeomorphism in- 
variance, one does not expect that such transformation 
will always respect the constraint We can overcome 
this problem, simply, by transforming the standard PPN 
metric to the gauge that satisfies ([24]) . This transforma- 
tion is given by 

5ij PN = 9ij ~ 2A 2 x,ij 

So7 N = 90i-(\i + X2)(Vi-Wi) 

9oo PN = 9oo ~ 2A 2 {U 2 + <t> w - $2) -2X 1 (A + B- $ x ) , 

(46) 

where the expressions fo r are functions of ten PPN 

parameters as given in [17]. Comparing eqs. (|46[) and 
(|43|) we can read off the values of the gauge parameters 
Ai and A2, as well as the PPN parameters as given in 
table H 

The parameters 7 and /3 are Eddington-Robertson- 
Schiff parameters used to describe the classical tests of 
theories of gravity; namely the deflection of light, time 
delay and perihelion shift. The parameter £ is non-zero 
in any theory of gravity that predicts preferred-location 
effects. Also, a\, oti and 013 measure whether or not the 
theory predicts post-Newtonian preferred-frame effects. 



2 We would like to thank Clifford M. Will for bringing this point 
to our attention. 



3 The reader can refer to Il9t for a recent and extended review of 
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parameter 


value 


effect 


limit 


7 - 1 


-3a 


time delay 


2.3 X 1U 






light deflection 


A w 1 A -4 

4 X 1U 


p — l 


8b 
32 a 


perihelion shift 


j X 11) 






Nordtvedt effect 


1.6 X 1U 


t 

? 


■<s„ 

8 a 


earth tides 


1 A - 3 
1U 


Q?l 


U 


orbital polarization 


1U 


Q2 


U 


orbital polarization 


4 X 1U 




8" 


orbital polarization 


A i n — 2U 
4 X 1U 


Ci 


fa 




2 x lO - ^ 


C2 


IB ° 


binary acceleration 


4 x f (T b 


C3 


—a 


Newtons 3rd law 


fO x 1(T S 


C4 









TABLE I: The values and limits on the PPN parameters 

When one attempts to devise integral conservation 
laws, we search for a quantity Q'*" which reduces to 
in flat spacetime and satisfies d^Q^ = 0. It was shown 
in [20( | that such quantity can exist only if all the pa- 
rameters {«3, Ci, C2, C3i Ci} are identically zero. Non-zero 
values of these parameters measure the extent at which 
a theory of gravity predicts violations of conservation of 
total energy and momentum. Notice that the parame- 
ter «3 plays a dual role, both as a conservation law and 
preferred-frame parameter. 

A bound ona 3 , of 4 x I0 -20 was reported in [2l|, |22j 
from the period derivatives of 21 millisecond pulsars. 
This small bound puts severe constraint on the value of 
{a<i}, and in turn on the diffcomorphism-violating La- 
grangians. 

Finally, it is worth noting that our result for goo to 
0(4) contains additional potentials UijUij and $5 that 
are not present in the standard PPN formalism. It would 
be interesting to devise experiments aimed to measure 
the effects of such terms in gravitational systems. 

10. SUMMARY 

Our work helps to quantify the physical content of gen- 
eral covariance. We have allowed for the possibility of 
small violations of diffeomorphism invariance through a 
class of operators with two derivatives of the metric. An 
analysis to linear order produced modest constraints from 
light bending. However, we have used the PPN formalism 
to bound the non-invariance of a sample operator in this 
basis, which produces a far stronger constraint. By far 
the strongest result comes from the absence of preferred 
frame effects in pulsars, and leads to the constraint that 
the dimensionless parameter a must be less than I0~ 20 
of gravitational strength. 

Tests of diffeomorphism invariance are of interest in 
its own right. We want to have quantitative probes of 



the different tests of GR. 



this fundamental property. Moreover, we have argued 
that this constraint is relevant for theories in which gen- 
eral relativity is an emergent phenomenon from a more 
fundamental theory that lacks a fundamental version of 
diffeomorphism invariance. There are many directions 
that extensions of this initial investigation can go and we 
feel that the topic deserves further study. 
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APPENDIX A: OPERATORS OF THE 
LINEARIZED THEORY 

In this appendix we write down the lower, marginal 
and higher dimensional operators for a linearized theory 
of gravity where we expand = rj^ v + . 

• Dimcnsion-2 operators 

The lowest dimension operators can only contain 
powers of the field, without derivatives. At leading 
order, 

h 2 ,h^h^ (Al) 

The Pauli-Fierz mass term is 

2 

— {h^h^-h 2 ) (A2) 

• Dimension-3 operators 



h\hhpl,KKhl (A3) 

• Dimension-4 operators 
At dimension-4 we start to have the possibility of 
operators with two derivatives, starting the next 
series in the derivative expansion. Those without 
derivatives are 



h\h 2 h^,hh^hl,h a ^h]hi (A4) 

while those with two derivatives are the set 

C( 4 ) = {d^d^ha/s , d^d-yh} , d a hd a h , dahd^h^} 

(A5) 
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• Dimcnsion-5 operators 
At this order we stop listing the series with zero 
derivatives and show the next order terms in the 
series with two derivatives. These all have three 
powers of the field and occur in the combinations 



C 4 h, C^ afi h aP 
where C 4 is defined above and 



(A6) 



d v h^d a h (3 \ 



(A7) 



• Dimension-6 operators 
At sixth order, the four derivative series starts. We 
here list only those with four derivatives and two 
powers of the field. 

d a d fi hd a d f3 h , dadphpudPdPhP" , d a d h^d t id u h af3 , 
d a dph llv &*ff , h% . (A8) 



These potentials satisfy the differential relations 
V 2 Vi = -Anpv, , Vi,i = -(7,o 

V 2 $l = -47T/W 2 , V 2 $ 2 = -4wpU 

V 2 $ 3 = -4wpU , V 2 $ 4 = -4?rp 
V 2 ($ w + 2U 2 - 3$ 2 ) = 2 X , ij U ij 
X,oo=A + B-* 1 , (B2) 

where 

X = - J d*x'p'\x-x!\ 
X,ij = SijU + Uij, V 2 X = -2U. (B3) 
In addition, consider the potential 



[ j3 / U ijP',j 

J \x-x'\' 



(B4) 



such that y 2 ipi = —AnUijpj. Hence, we define the po- 
tential $5 as 

V 2 i/> M = V 4 $ 5 = 4irp ti U ti - AnUijPM . (B5) 



APPENDIX C: EXPRESSIONS USED 
THROUGHOUT THE PAPER 



APPENDIX B: GRAVITATIONAL POTENTIALS 



In this appendix, we list the various gravitational po- 
tentials used to construct the metric. 



U, 



J \x-x'\ 

,P' (x-x')i [x-x') j 



w 



d A x' 



\x — x 



•'13 



J \x-x'\ 

J \x-x'\ J \x-x' 

$3 = [d 3 x'-^-,^= fd^x'-P^ 

J \x — x'\ J \x — x' 

I 



B = d 



, 3 , p' 'dv 1 1 dt ■ (x — x') 



\x — x' 



d A x'p'p" 



X ' — x' I 3 \\x — x"\ \x'—x 



'ii \™i _ ■ 



(Bl) 



In this appendix we give the form of the different ex- 
pressions used throughout this paper. 

The components of the affine connections are 



(2) 

V 1 - 

1 oo — 

(3) ■, 

r l — - 

L 0j - 9 



1 (2) 



(4) 



1 (4) , (3) 1 (2) (2) 



ff00,i , r 00 — — - g 00}i + 5oi,0 + " 9ij900,j 



(2) 

(4) 
-pi 

L jk ~ 



(3) (2) (3) 

9i0,j + 9ij,0 — 9j0,i 



(2) (2) _ (2) 

9ij,k T" 9ik,j 9jk,i 



(3) 



2 

1 (2) 



^oo— 2^o°>° 

(2) 1 (2) 

L 0i— ~2 ^ 00 ' i 



-S ip 
2 



(4) (4) _ (4) 

9pj,k i 9pk.i 9ik,p 



l (2) 
9ip 



(2) (2) _ (2) 

lpi,k T" 9pk,i 9ik,p 



(Cl) 



The lower order expressions for R and M. read 



(2) 

-Roo 

(2) 

Rij 

(2) 

Mi, = - 

(3) 

Rio 



1 (2) (2) 3 (2) 

— 2 900,ii , -Moo— ^ 900, i: 



(2) _ (2) (2) (2) _ (2) 

900,ij 9kk.ij T 9ik,kj T" 9kj,ki 9"ij-,kk 



(2) _ (2) (2) 
9ik,kj 9jk,ki TJ 9ij,kk 



(2) (3) (2) (3) 



(3) 



3 (2) 

2 5i0,fcfc 



(C2) 



10 



The higher order perturbations of R and A4 are given 
by (the constraint (|36[) being imposed) 



(4) 

Roo = 



(4) 

Moo 



(4) 



1 (2) (3) 1 (4) 1 (2) (2) 

— 2 9m, 00 + SiO.iO ~2 300,fcfc +^ 9ij900,ij 



1 (2) (2) 
T <?00,i<?00,i 



1 (2) (2) 
•7 900,i9jj,i 



15 (2) (2) (2) (2) (2) (4) 

9oo,i9oo,i +3 gujgooj —2 300,00 +2 ffoojj 

(3) (2) (2) 3 (2) (2) (2) (2) 

-4 goj.oj — gijgoo.ij gij,kgij,k — gjk,igij,k 

3 (2) (2) 3 (2) (2) 3 (2) (3) 
"J 9ii,j900,j +^ 9kk,j9ii,j ~ 2 fl^OO + <70i,0i 

3 (2) (2) 9 (2) (2) 3 (2) (2) 
"g #00,i<?00,i +g gik,jgik,j ~ J 9ij,k9kj,i 

3 (2) (2) 1 / (4) (4) \ 

"2 9kp9ii,kp +2 I 5ii,fcfe <?ifc,ifc I ) (^oj 



(4) 

where the summation is implied in .M^ 
The functions and Pjj defined in 



the trace of {7| , and write h in terms of T to find 

(1 + 01 - 3a 3 ) □/i Q/3 + (1 - a 2 + 2a 4 ) d a d h 

( (2) (2)\ 
= -WnG T af} -Ari afi T 



where 



A = 



1 — 02 + a 4 + a 5 



2 — ai — 3a2 + 3a 3 + 2a 4 + 4a5 
Solving for the hoo component we get 
hoo = 2aU , 

where 

2G (1 — a! — 2ci2 + 3a3 + a 4 + 3as) 



(Dl) 



(D2) 



(D3) 



a 



(1 + 01 — 303) (2 — ai — 3a2 + 3a3 + 2a 4 + 4as) 



are given by 



(D4) 

The most general solution of is given by 

hij = axSijU + (72?/^ , (D5) 



(4) 

Bh = 



(2) 



(4) 

Q fe 



- ,'3 (2) (2) 3 (2) (2) 
— I o 300,fc900,fe +0 <?nfc,m<?nfc,; 



1 (2) (2) 
- 9 nrn, k 9 km 7 7i 



(2) (2) 



(2) (2) 



1 00 1 00+1 fern 1 fem 



+2 r 



(2) (2) (2) (2) 

m-pm 1 o -pO "pO 
ik l jk ~r z 1 iO 1 jO 



(2) (2) 



(2) (2) 



(3) 



r° A k 4- r m A k 4- 4° 



(4) 



1 (2) 
" 2 9ip 



(2) 



(2) 
9pk,j 



(2) 
9jk,p 



1 (2) / (2) (2) _ (2) 

2 I 9pi,k "T" gpk,i gik,p 



1 (2) 



(2) 
9pi,j 



(2) 



(2) 



(2) 



(2) 
9pj,i 

(2) 



(2) 

9ij;P 



(2) 



(2) 



+ 0niP?fc + ffnjltf. - 3fe m r 



W 1 im 



(C4) 



where cti and 02 are constants to be determined. Using 
the constraint dihij we find a\ = a-i- Then, substituting 
into eq. (|D1[) we find that the following two equations 

4G(l-a 2 + a 4 + a 5 ) 



2<Ti(l +a± — 3a 3 ) 



and 



2 — ai — 3a2 + 3a3 + 2a 4 + 4as 

(D6) 



<7 1 (—1 + ai — 3a 3 + 2a 2 — 4a 4 ) = —a (1 — a 2 + 2a 4 ) , 

have to be satisfied simultaneously. This can be true only 
if we take a 4 = 05. This is exactly what we found before 
in eq. ©. Setting a — 1, we obtain the normalized value 
of a 1 



0\ 



1 — Gi2 + 2a 4 



(D8) 



1 — a% — 2a2 + 3a3 + 4a 4 
from which we immediately read the PPN parameter 7 

1 — (22 + 2<2 4 
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1 — Oi — 2a2 + 3ci3 + 4a 4 



(D9) 



APPENDIX D: LINEARIZED VERSION OF THE 
EQUATIONS OF MOTION 



which reduces in the limiting case a\ = 0,2- ■■ = 05 = to 
the GR result. The deflection of light is proportional to 
7 + 1, which gives (fTT)) . 



One can use the linearized version of the equations of 
motion (J7J) along with the constraint dph a P — to solve 
for gij to 0(2) and goi to 0(3). 4 To this end, we take 



4 Notice that hoo = <?qq , hij 



■ g\f , and h oi 



(3) 



APPENDIX E: SOLVING FOR I 3 g-tlk -2 glt]ik 



To find g^ to 0(4), we need to write the i — j com- 
ponent of the equation of motion ([33]) to 0(4). How- 



ever, since we are interested in 

-,(4) 



(4) 



9oo 



to first order in 



a, we can solve for g\- to zeroth order in a. Hence, 
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the i — j component of (|33| reduces to the GR result 
Rij+0{a) = 8irG {T l3 - g l} T/2). Moreover, the solution 

of g!-j should respect the constraint ([55]) . 

The i — j component of the Ricci tensor can be evalu- 
ated to the fourth order to find 

(4) (4) (4) (4) 

gij,kk + 9kk,ij — gik,kj ~ 9jk,ki— ^ij l m J 

where Sij is a complicated function of the various poten- 
tials. Contracting eq. (|E1[) results in 



(4) 

The value of 500 can be calculated using R 00 + 0(a) = 
8irG (Too — gooT/2), and imposing the constraint (136)) to 
find 



500= -U 2 + $w + 4$i - 4$ 2 + 2$ 3 + 6$ 4 ■ (E8) 



Hence, using the dictionary in appendix F we finally ob- 
tain 



g\ 



(4) (4) 



1 



(E2) 



while differentiating it with respect to j gives 

(4) (4) \ _ , . 

9ij,ij I — ^mk.k ■ {^"^J 



V 2 g\ 



From eqs. (|E2[) and (|E3|) we obtain the intcgrability con- 
dition 



25 5 15 

— U 2 + ~®w + 8$l + y$2 + 8$ 3 



1, 



(E9) 



Writing i£j as the sum of the gradient and curl of a scalar 
and vector, i.e. E.- L = V,, + (V x 4)„ we find 



Sij ^SijSkk 



0. 



(E4) 



(E10) 



It was shown by Chandrasekhar and Nutku [T|| that this 
condition is indeed satisfied in GR in the PPN gauge. 
Since GR is a diffcomorphism-invariant theory, we con- 
clude that this condition still holds when using the con- 
straint (|36f . It was also shown in [] that the solution of 
eq.(|E2|) is given by 



(E5) 



V 2 gij= S i3 + E it j + Ej^ , 



where Ei are arbitrary functions. 5 Now, using eqs. (|E2 
and (|E5[) we get 



APPENDIX F: DICTIONARY 



In this appendix, we give a dictionary for the different 
combinations that appear in our formalism. 

Using the differential relations in appendix B, we ob- 
tain to the zcroth order of a 



(4) (4) 

3 ga^kk 2 g%k^%k— 2Sa + 2Eij 



(E6) 



where 



1 (2) (2) (2) (2) 3 (2) (2) (2) (2) 

bii — 2 S00,i<?0CM + 900900,ii + Z 9kp,i9kp,i + fffcpfffcp,zi 

(2) (2) (2) (2) (2) 1 (2) (2) 

+ ^ 9ii,00 + 9kp9ii,pk ~ 9ik 



, 1 (2) (2) (3) (4) 

+ 2 900,k9ii,k goifii + .900,fcfc 



(4) 



(4) 



(2) 



(2) 



I* +3 T 00 -3 St 00 + St 00 1 . (E7) 



5 Strictly speaking, since we are using only GR, the constraint (136 I t 
is not more than a gauge fixing to 0(2), and Ej accounts for the 
gauge freedom to 0(4) . Although we are free to fix a gauge and 
hence the values of Ei in case we were dealing only with GR, 
the values of these functions will be determined upon using the 
consistency condition i !24[ ) to 0(4). 



(2) (2) 
Qij,kQij,k 

(2) (2) 
9ij9ij,kk 

(2) (2) 
9jk9ii,jk 

(3) 
90iM 

(2) (2) 

3oo, i.9oo, i 

(2) (2) 
9ii 7 k9jj?k 

(2) 

3oo,oo 

(2) (2) 
9ik,j 9ij,k 

(2) (2) 
9oa,i9kk,i 

(2) (2) 
.900.900, kk 

(2) (2) 
9ijg00,ij 



= V 



~U 2 



3> 



w 



V 2 [-2C/ 2 - $ w - 
V 2 [4C/ 2 + 2<±V H 
V 2 [$i-A- B] 
V 2 [2C/ 2 - 4$ 2 ] 
V 2 [8C/ 2 - 16$ 2 ] 

V 2 [*! - .A - 6] 



= V 



;c/ 2 



7*3 + ^«C^ 

f 7$ 2 ] 
-2* 2 1 



$2 + ■zUuUi. 



V 2 [4J7 2 - 8$ 2 ] 
V 2 [4$ 2 ] 

V 2 [2U 2 + $w + $2] 



(Fl) 
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One can also show 

U,<jU,ij = V 4 (hj 2 - i$ 2 j +iTTp, i U, i 

/3 1 5 
X,ijkU,ij k = V 4 ( -t/ 2 + -$w - -$ 2 



X,ijkmX,ijkm = V 4 ( — [/ 2 + $w - 4$ 2 + ^E/yE/y 



+4wp ii U t i + 2nU lj p, i: j - 2nUV^p 

'11 9 1 
— (7 2 + $w - 4$ 2 + t* 

+8Trp,iU ti + 4arUijp tij - 47rJ7V 2 p 
V 4 $ 2 = -8irp,iUi + 16n 2 p 2 - 4nUV 2 p 

V 4 $ 4 = - 1/3] - inpjUi + 16ttV . 



(F2) 
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